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Abstract 

We investigate theoretically the photoacoustic generation by a gold nanosphere in water in 
the thermoelastic regime. Specifically, we consider the long-pulse illumination regime, in which 
the time for electron-phonon thermalisation can be neglected and photoacoustic wave generation 
arises solely from the thermo-elastic stress caused by the temperature increase of the nanosphere 
or its liquid environment. Photoacoustic signals are predicted based on the successive resolution of 
a thermal diffusion problem and a thermoelastic problem, taking into account the finite size of the 
gold nanosphere, thermoelastic and elastic properties of both water and gold, and the temperature- 
dependence of the thermal expansion coefficient of water. For sufficiently high illumination fluences, 
this temperature dependence yields a nonlinear relationship between the photoacoustic amplitude 
and the fluence. For nanosecond pulses in the linear regime, we show that more than 90 % 
of the emitted photoacoustic energy is generated in water, and the thickness of the generating 
layer around the particle scales close to the square root of the pulse duration. The amplitude 
of the photoacoustic wave in the linear regime are accurately predicted by the point-absorber 
model introduced by Calasso et al. [], but our results demonstrate that this model significantly 
overestimates the amplitude of photoacoustic waves in the nonlinear regime. We therefore provide 
quantitative estimates of a critical energy, dehned as the absorbed energy required such that the 
nonlinear contribution is equal to that of the linear contribution. Our results suggest that the 
critical energy scales as the volume of water over which heat diffuses during the illumination 
pulse. Moreover, thermal nonlinearity is shown to be expected only for sufficiently high ultrasound 
frequency. Finally, we show that the relationship between the photoacoustic amplitude and the 
equilibrium temperature at sufficiently high fluence reflects the thermal diffusion at the nanoscale 
around the gold nanosphere. 
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I. INTRODUCTION 

Photoacoustic imaging is a promising modality for biomedical applications that has 
emerged during the last two decadesii^. This non-invasive modality is based on the con¬ 
version of absorbed light energy into ultrasound via the thermoelastic effect^. The image 
contrast therefore depends on the optical absorption properties of the medium. In biological 
tissues, endogenous optical absorption may be used to form various types of images. For 
instance imaging of the hemoglobin enables reconstruction of the vascularization network-. 
To further enhance the contrast or obtain complementary information, various exogenous 
contrast agents have been developed for photoacoustics^. In particular, plasmonic noble 
metal nanoparticles have been introduced as photoacoustic contrast agents in the early 
20008^“^. Gold nanoparticles (GNP) are very attractive as photoacoustic contrast agents 
thanks to their large optical absorption cross-section, their resistance to high illumination 
fluences and their spectral selectivity based on their surface plasmon resonance^. Typically, 
the optical absorption cross-section of noble metal nanoparticles is a few orders of magni¬ 
tude larger than that of traditional molecular dyes^*^. One major consequence of the strong 
optical cross-section of GNPs is that their temperature can increase signihcantly when they 
absorbe pulsed light^, leading to various possible phenomena including nano/micro-bubble 
formatiort^^— or even modifications of the particle shape. In the context of photoacoustic 
imaging, bubble formation is interesting as the emitted signals are usually strong and exhibit 
a nonlinear relationship between the incident fluence and the photoacoustic amplitudei^i^. 
More generally, nonlinear photoacoustic phenomena provide a means of selectively detecting 
contrast agents from an absorption background that behaves linearly!^, similarly to what 
is done in the held of ultrasound imaging. Several phenomena may induce nonlinear rela¬ 
tionships between the photoacoustic signal amplitude and the energy of the incident light 
in addition to bubble formation, such as optical saturation^ii^, photochemical reaction^^ 
and temperature-dependent thermodynamic parameters^'i^'i^^— . The latter phenomenon 
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is related to the temperature dependence of the thermal expansion coefficient, and is at 
the core of this paper, in which we investigate theoretically the photoacoustic generation 
by a gold nanosphere in the thermoelastic regime (absence of bubble formation). In this 
work, our study is restricted to what we call the long-pulse illumination regime (typically 
longer than several picoseconds), for which the typical electron-phonon relaxation time in 
the gold nanoparticle is negligible compared to the illumination duration. In particular, 
electrons and phonon may be considered not only to have thermalized distributions (such 
as assumed in the two-temperature model^), but also to be described by a unique tem¬ 
perature value, a situation very different from that encountered with sub-picosecond light 
pulses where nonthermal distributions may be involved2^>^. In the long-pulse illumination 
regime, a single value for temperature may be used to describe thermodynamic properties, 
and the photoacoustic wave generation arises from the thermo-elastic stress caused by the 
temperature increase of the nanosphere or its liquid environment. This regime encompasses 
the nanosecond-pulse regime commonly used in photoacoustic imaging, and to which most 
of our study of thermal nonlinearity is restricted. 

Although nanoparticles have now been used for more than a decade as contrast agents for 
photoacoustic imaging in many research studies, comparatively few studies investigated the 
physics involved at the scale of the nanoparticle, whether with theoretical or experimental 
approaches. 

In the linear regime, Inkov et ai— introduced a three-step model to predict the pho¬ 
toacoustic emission by an absorbing spherical particle, based on solving (1) a light absorp¬ 
tion problem, (2) a thermal diffusion problem and (3) an acoustic problem. Analytical 
solutions were provided in the linear regime, but were limited to thermally small or large 
particles, assumptions that are not valid for gold nanospheres illuminated with nanosecond 
pulses, as commonly encountered in photoacoustic imaging. For nanospheres illuminated 
with nanosecond pulses, the thermal relaxation time is comparable to the pulse duration, 
making it complex to analytically derive the temperature held. Numerical approaches are 
usually required to solve the thermal problem in this case, as was done for instance by Baf- 
fou and Rigneault— . More recently, it has been demonstrated experimentally in the linear 
regime that for gold nanospheres illuminated with nanosecond pulses, it is mostly the liquid 
surrounding the nanoparticles that emits photoacoustic waves^I, as was theoretically dis¬ 
cussed earlier for thermally small particles^. A theoretical explanation in the linear regime 
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was proposed by Chen et al.—, but the temperature held was modeled via a quasi-static 
thermal held. Other studies provided theoretical expressions for the photoacoustic emission 
by spherical absorbers, but with assumptions not always valid in this work and more impor¬ 
tantly also limited to the linear regime2^>2£. Two recent studies also reported comparisons 
between experimental results and theoretical predictions in the linear regime^ 1 ^, but the 
authors assumed heat and stress conhnement at the scale of a whole particle suspension and 
did therefore not consider the photoacoustic generation at the scale of individual particles. 

In the nonlinear thermoelastic regime, a few works have been reported with nanoparticles. 
Calasso et al. - provided analytical expressions of the waveforms emitted by a point-absorber 
model. In our work, we extensively use the predictions of this model to be compared to 
our own model that takes into account the hnite size of the gold nanosphere, both in the 
linear and nonlinear regime. Following the physical approach introduced by Inkov et al.—, 
Egerev and Oraevsky— described and observed the nonlinear photoacoustic generation by 
gold nanospheres in water in the thermoelastic regime. A simple criteria to assess the 
signihcance of nonlinear generation was derived based on scaling arguments, that shows that 
thermal nonlinearity should be observed at sufficiently high fluences (Eq. (3) of Egerev and 
Oraevsky—), as demonstrated experimentally. The validity of this criteria will be discussed 
in sections IIII B II and IIIIB41 in comparison with the quantitative predictions from both 
the point-absorber model and our model. In the thermoelastic regime, another previous 
experimental study reported a nonlinear increase of the photoacoustic signal with the laser 
fluence^. The origin of the nonlinearity was assumed to be the temperature-dependence 
of the thermal expansion coefficient, but presumably caused by thermal coupling within 
aggregated nanoparticles in cells, rather than by temperature elevation around individual 
nanoparticles as considered in our work. 

The main objective of our theoretical work is to provide physical insight and quantitative 
predictions regarding the photoacoustic generation by a single gold nanosphere, for both 
the linear and nonlinear thermoelastic regimes, beyond the few initial results from previ¬ 
ous relevant studies^'^'^'^. In section m we hrst introduce the physical model used to 
predict the photoacoustic signal from a gold nanosphere. Our model takes into account 
the hnite size of the nanosphere and its elastic and thermo-elastic properties, as well as the 
temperature-dependence of the thermal expansion coefficient of the surrounding liquid. Pre¬ 
dictions from the point-absorber model are also given for further comparison in the results 


6 







section. Section mi also describes the principles that were used for solving the thermal and 
acoustic problems, with details further given in Appendixes |B] and O Section [TTIl provides 
results and discussions in both the linear and nonlinear thermoelastic regimes. In the linear 
regime, we study the origin of the photoacoustic wave as a function of the size of the gold 
nanosphere and the pulse duration. From there, the paper focuses on the nanosecond pulse 
regime for which the generation mostly occurs in the surrounding liquid. A scaling law is 
found that describes the typical thickness of the water layer that generates the photoacoustic 
wave. A comparison between our results and those from the point-absorber modeU^ is then 
provided, as a preamble to further analysis in the nonlinear regime. In the nonlinear regime, 
we hrst derive predictions from the expressions given by the point-absorber model. We 
then describe quantitatively temperature rises in gold nanospheres, showing that thermal 
nonlinearity is indeed expected at commonly encountered light fluences. Our predictions 
for the photoacoustic amplitude emitted by a gold nanosphere in the nonlinear regime are 
given and compared to those from the point-absorber model, and we study the occurence of 
nonlinearity as a function of fluence and particle size. Finally, we investigate the influence of 
the equilibriume temperature on the photoacoustic amplitude, for a given nanosphere and 
hxed fluence. 

II. THEORETICAL APPROACH 

We consider the photoacoustic generation by a single gold nanosphere in water. In this 
section, we describe the models and computations that are used to produce the theoretical 
predictions given in Sec. IIIII Throughout all this work, spherical symmetry is assumed, with 
the center of the gold nanosphere as the center of symmetry. 

A. Physical models 

1. Photoacoustic generation by a point absorber in a liquid 

In order to discuss further below the generation of photoacoustic waves by a gold 
nanosphere and to introduce relevant quantities used throughout our work, we hrst re¬ 
call here the basic equations which describe the photoacoustic generation in the simple 
situation of a mechanically and thermally homogeneous liquid medium. In particular, we 
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provide the analytical expressions derived by Calasso et ai— of the photoacoustic pressure 
wave emitted for the limiting case of a point-absorber, both in the linear and nonlinear 
thermo-elastic regime. The corresponding analytical predictions are discussed in the results 
section, in particular in comparison with the theoretical predictions for a gold nanosphere. 
We will show that the point-absorber model accurately predicts the photoacoustic emission 
by a gold nanosphere only in the linear regime, which justihes the introduction of our 
nonlinear model for the gold nanosphere in the next section. 

a. Model equations. When the physical properties are assumed to be homogeneous and 
constant in time, the generation of photoacoustic waves in a liquid medium is dictated by 
the following system of coupled equations^'i^'^^ 

dT 

PoCp-^{r,t) - KAT{r,t) = Py{r,t) (1) 

1 d‘^v d'^T 

Ap(r, t) - ^^(r, t) = -po/3o^(r, t) (2) 

where p{r,t) and T(r, t) are respectively the (photo)acoustic pressure held and the tem¬ 
perature held. The relevant physical properties are the mass density po, the coefficient of 
thermal expansion (3o = —^ (ot)o’ acoustic velocity Cg, the thermal conductivity k and 

the specihc heat capacity at constant pressure Cp. Eq.(IT]) is a standard heat dihusion equa¬ 
tion, with a heat source term P^(r, t) representing the volumetric density of power converted 
to heat (with dimensions unit power per unit volume). Eq.([T]), as fully decoupled from the 
pressure held, is only valid for liquid and solid media (as opposed to gases) for which the 
ratio 7 = — of the specihc heat capacity at constant pressure to the specihc heat capacity 
at constant volume may be considered very close to 1 (see demonstration in Appendix lA II 
based on Morse and Ingard— ), an assumption that will be made throughout all this work. 
The photoacoustic wave equation ([2]) is a classical wave equation with a source term given 
by the second time-derivative of the temperature held. 

Equations ([T]) and (EJ indicate that to solve the photoacoustic problem given a source 
term Py{r,t), one has to first solve a thermal dihusion problem, and then solve an acoustic 
problem once the source term given by the temperature held is known. In the context 
of photoacoustics, the heat source term arises from optical absorption, and is therefore 
proportional to some illumination function such as the huence rate (or intensity) <hr(r,t) 
(unit power per unit surface). In this work, we will considered a single optical absorber 
(with an absorption cross-section cXa) illuminated with some incident pulsed light described 






by the following expression 


( 3 ) 


'^P '^P 

where $o is the fluence (unit energy per unit surface) and / is a dimensionless peaked 
function describing the temporal prohle of the fluence rate. / verihes /j!(^/('f)df = 1 and 
is normalized such as Tp is dehned as the full width at half maximum (xp is further referred 
to in the text as the pulse duration). Throughout this work, the laser temporal prohle is 
chosen as a Gaussian dehned accordingly by 


/(^) = 


2-\/ln(2) ^_ 4 in(- 2 )T 


TT 


( 4 ) 


b. Point absorber in the linear regime. For an optical absorber of vanishingly small size, 
but with a hnite optical absorption cross-section (Tq, Calasso et al.'^~ provided an analytical 
expression of the photoacoustic pressure wave emitted by the ’’point-absorber” (referred to 
as a photoacoustic point source in Calasso et al.—)\ 

1 1 d/ 


p(r, t) = EahsPo 


CpTp dvrr dr 


T = 


t-^ 

Cs 


( 5 ) 


where Eats = o'a'ho is the energy absorbed by the point absorber. As expected from the 
linearity of Eqs. ([I]) and ((2]), the photoacoustic pressure is proportional to the absorbed 
optical energy. 
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FIG. 1: Thermal expansion coefficient of water as a function of temperature. Data derived 
from the density of water as a function of temperature^ 
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c. Nonlinear thermo-elastic regime. When significant temperature rises occur, the 
physical properties involved in ([1]) and (|2]) may vary during the illumination and subsequent 
photoacoustic generation. It is well known that amongst the relevant thermodynamics 
properties, the thermal expansion coefficient [3 shows the most signihcant temperature de¬ 
pendency^. The temperature dependence of the thermal expansion coefficient /3(T) of water 
is shown in Fig. [1] Taking into account this temperature dependence, Eq. [2] has to be 
modihed as 

c«Ap(r,t) - = -Pocl^ 

Whereas the temperature held remains linearly related to the optical illumination (via Eq. [1]), 
the photoacoustic pressure wave in Eq. [6] is nonlinearly dependent on the temperature held 
and therefore on the optical illumination. Note that the source term in Eq. [6] is slightly 
diherent from that given initially in the pioneer work by Burmistrova et ai— or in Gusev 
and Karabutov- (see detailed derivation of our equation in Appendix lA 1|) . 

d. Point absorber in the nonlinear regime. For the point absorber model, under the 
assumption that the temperature dependence of the thermal expansion coefficient can be 
linearized as /3{T) = Po + /5i(T — Tq) (with Po = /3{To) and (3i = ^(Tq)), and for a 
gaussian temporal prohle of the illumination pulse, Calasso et ai— also provided an analytic 
expression of the photoacoustic pressure wave emitted by the point absorber in the nonlinear 
regime: 


p(r, t) =EabsPo 


1 1 d/ 

CpTp dvrr df 



EL/3i- 


, 3 / 2 ^ 2 , 


- 7/2 4 ^^ 


h 



(7) 


where h(r) is a dimensionless function with a tripolar shape, given in details by Eq. flA16D in 
Appendix IA3I Note that the numerical prefactor in the nonlinear term of Eq. [7] is diherent 
from that the original equation (25) given in Calasso et al.—, for various reasons detailed in 
Appendix lA 11 including the modihcation required to take into account the correct source 
term of Eq. |6l 
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2. Photoacoustic generation by a gold nanosphere in a liquid 


In this section, we present the physical models used to describe the photoacoustic emission 


by a hnite-size gold nanosphere (of radius Rg) immersed in water. As opposed to the case of 


a point absorber, no analytical expression is available for the emitted photoacoustic pressure 
wave by a hnite-size solid sphere, except for very limiting cases in the linear regime^^i^Si^^ 
out of scope here. In Chen et al.—, the thermal source term was modelled via a quasi¬ 
static thermal held, and the solution in the linear regime was approximated in the Fourier 
domain assuming the sphere was small compared to the ultrasound wavelength. The model 
used in our work takes into account thermal dihusion around the nanosphere, photoacoustic 
generation and propagation in both the gold nanosphere and its liquid environment, and 
any arbitrary temperature-dependence of the thermal expansion coefficient of the liquid 
environment. 

a. Thermal model. Because the thermal conductivity of gold is much larger than that 
of water, the temperature within the nanosphere is considered in this work to be uniform, 
which is known to be a very accurate approximation for gold spheres of diameter of the 
order of a few tens of nanometers, and for pulse duration no shorter than a few In 

particular, we emphasize that for the pulse durations considered in our work, longer than 
a few tens of ps, the electron-phonon thermalization that occurs on a time scale no longer 
than a few ps can be totally neglected^i^^, as opposed to the situation usually encountered 
in picosecond acoustics with sub-picosecond illumination^^*^^. Under this assumption, the 
spatio-temporal evolution of the temperature helds Ts{t) inside the solid gold nanosphere 
and Tf{r, t) in its liquid environment can be described by the following system of differential 
equations^^: 



PoCp 


(8b) 


(8a) 


with the following boundary conditions: 


n(fi.+,i) = r,(i) 

Tfl'r co.t) ^ To 


(9a) 

(9b) 
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with r = ||r||, and the subscripts s and / referring respectively to the solid and fluid 
phases. Eq. [Ha] states that the variation of the uniform sphere temperature increases via the 
absorbed optical energy and decreases via thermal conduction at the gold/water interface. 
This equation takes into account the continuity condition for the thermal flux across the 
interface. Eqs. [Hal and |9b] provide the additional boundary conditions required to solve the 
problem. The continuity equation Eq. [9a| assumes that any interfacial thermal resistivity is 
neglected. This assumption is discussed further in section IIII Cl 

b. Thermoelastic equations. Under spherical symmetry and for isotropic materials, 
the thermoelastic equations in both heterogeneous solid and liquid media can be written as 
a hrst order velocity-stress system of equations that reads^^^^ 


dvr 

dt 




darr 




daee 

dt 




1 r darr / \ 2 

H-rw -^[r,t) +-[arr-(Tee) 

Po[r) dr r 


(A(r) 2p(r))^ 2A(r)- 

dr r 


Vr{Y,t) 


- (^(0 + \T{r))(3{T{Y,t))^{Y,t) 


A(r)|- + 2(A(r)+^(r))l 

dr r 


Vr{Y,t) 


2 dT 

-{X{r) + -p{r))/3{T{Y,t))^{Y,t) 


(10a) 

(10b) 

(10c) 


where a is the stress tensor, Vr is the radial displacement velocity, and A and /i are the Lame 
coefficients. One can readily verify that if p is set to zero in Eqs. HUl i.e. the material is a 
liquid (arr = —p), the system yields Eq. [6] (or Eq. |2]for constant (3). In the relevant case 
here of a solid/liquid interface, the following continuity conditions must hold for both the 
radial velocity and stress at the sphere interface : 


Vr{Rs ,t) = Vr{Rt ,t) (11a) 

arr{R~R) = -p{Rt, t) (11b) 


B. Computations for a gold nanosphere 

The equations that describe the photoacoustic generation by a solid and optically ab¬ 
sorbing sphere (Eqs. IHlto fTTl) are much more complex than the equations for a homogeneous 
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liquid (Eqs. [T] to [2]) and cannot be solved analytically. However, their resolution still re¬ 
quires to first compute the temperature held from the thermal problem, and then to use 
this temperature held as a source term in the thermoelastic problem. The full resolution of 
both the thermal and thermoelastic problems is referred to further in the text as a numerical 
simulation, based on the computational approaches described below. 


1. Temperature computations 


The system of equations [H] and [H] may be solved analytically for an impulse excitation, 
i.e —)■ S(t), using the Laplace Transform. After tedious but simple algebric ma¬ 

nipulations and use of tables of known inverse Laplace transforms, one may obtain the 
Green’s function Gth(r, f) (solution to a 6{t) excitation) of the thermal problem, as was 
done by Egerev et al.—. The expression of Gth(r, t) is given in Appendix [Bl The temper¬ 
ature held in water for a pulse excitation can then be calculated by the convolution of the 
thermal Green’s function with the source term: 


r(r,t) = G'th(r,f) * 


O'abs^O 




Tp Tp 


( 12 ) 


For all our results, the convolution in Eq. [T2] was performed numerically, with the source 
function / given by Eq. 01 The temperature held T{r,t) was computed and sampled on a 
regular grid T{n x Ar,m x At) required by the hnite-diherence in time-domain resolution 
of the thermoelastic problem described below. 


2. Acoustic computations 

In this work, we used a hnite-diherence time-domain (FDTD) algorithm to solve the 
thermo-elastic problem. We adapted the well-known Virieux’s scheme to our problem with 
spherical symmetry. In brief, the Virieux’s scheme for elastodynamics^ (analog to the Yee’s 
scheme for electromagnetisnt^^) is based on a spatio-temporal discretization of the system of 
continuous equations (Eqs. [TD]) on staggered grids (spatial grid step Ar and temporal grid 
step At). The solution is computed step by step in time, over the whole spatial domain at 
each time step. Any known source term may be taken into account, both in the sphere and 
in water. In particular, it makes it straightforward to take into account the temperature- 
dependence of the thermal expansion coefficient of water, by simply computing the value 
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TABLE I; Physical constants associated with gold and water, at T 25°C, from^ 


Properties Gold Water Unit 


Mass density po 
Specific heat capacity Cp 
Thermal conductivity k 
Thermal diffusivity y 
Thermal expansion /3 
First Lame coefficient A 
Second Lame coefficient p 
Compressional wave velocity 
Shear wave velocity 


19.3 

1.00 

xlO^kg.m ^ 

129 

4200 

J.kg-hK-i 

318 

0.60 

W.m-hK-i 

128 

0.142 

X 10“®m^.s“^ 

0.43 

FigU 

xlO-^K-i 

147 

2.25 

GPa 

27.8 

2.25 

GPa 

3.24 

1.50 

pm.ns“^ 

1.20 

- 

pm.ns“^ 


of P^{Tf{r,t)) at each point in space and time. Another well-known key advantage of 
the Virienx’s scheme is that bonndary conditions snch as given by Eqs. [H] are implicitly 
verihed^. As a conseqnence, reflected and transmitted aconstic waves at the water-gold 
interface were taken into acconnt in our numerical solutions. The discretized equations that 
were used are detailed in Appendix O The spatial grid step Ar was chosen small enough 
to ensure a proper convergence of the FDTD solution: the convergence was ensured by 
verifying that results with two different spatial steps showed no signihcant difference. The 
values of Ar typically ranged from 0.1 nm to 5 nm depending on the sphere radius and the 
pulse duration. The dimension of the spatial domain was taken sufficiently large (typically 
several tens of pm) such that any spurious reflections from the domain boundary would 
arrive far after the photoacoustic pressure waveforms. The time step At was derived from 
Ar via the stability condition given in Appendix O 


3. Values of the physical properties 

All the values of the physical properties of gold and water used in the computations are 
summarized in Table [T1 Except for the thermal expansion coefficient whose value may depend 
on temperature, the values for all other properties (assumed to be constant) were those at 
room temperature (~ 25°C). The absorption cross-section of a gold nanosphere depends 
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on its size, and therefore so does the absorbed energy for a given fluence. The values of 
absorption cross-sections used in this work were derived from the Mie theory with optical 
constants from Johnson and Christy—, for an illumination wavelength A = 532 nm. A few 
typical values are given in Table HIl For a diameter below typically 50 nm, the absorption 
cross-section is much larger than the scattering cross-section and scales as the nanoparticle 
volume^. 

TABLE II: Absorption cross-section aabs of a gold nanosphere as a function of radius, for 

an illumination wavelength A = 532 nm 

(nm) 10 20 30 40 50 60 

aabs (xlO-1^ m2) 3.7 x 10-2 3.3 x lO’^ 1.1 1.8 2.1 2.3 


III. RESULTS AND DISCUSSION 

A. Linear regime 

In this section, the temperature-dependence of the thermal expansion coefficient of water 
is neglected, i.e we consider the linear photoacoustic regime. Our hrst objective is to in¬ 
vestigate the relative contribution of the gold nanosphere and its liquid environment to the 
sound generation. When the photoacoustic wave is predominantly generated from the liquid 
environment rather than from the solid sphere, we then investigate the typical thickness of 
the water layer that generates the photoacoustic wave. As a preamble to our results in the 
nonlinear regime, our results for the gold nanosphere are compared to those predicted by 
the point absorber model in the linear regime hrst. Throughout all the paper, all absolute 
photoacoustic amplitudes are given at 1 mm from the center of the absorber (r = 1 mm). 
In the linear regime, all the results are predicted for an equilibrium temperature Tq = 20°C'. 


1. Origin of the photoacoustic wave 

The absorption of the laser pulse by a spherical gold nanoparticle creates a transient 
temperature rise in both the particle and its liquid environment due to heat diffusion. 
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X = 15 ps X = 500 ps X = 5 ns 

P P P 



Time (ns) 


Time (ns) 


Time (ns) 


FIG. 2: Photoacoustic waveforms obtained in the linear regime, for a 20-nm radius sphere, 
illuminated with a fluence $ = 1 mJ.cm“^. The red curves correspond to signals obtained 
with thermal expansion in both water and gold (SI). The blue curves correspond to signals 
obtained with no thermal expansion in gold (S2). 


From Eqs. [101 if is clear that both the gold nanosphere and its environment may generate 
photoacoustic waves. Here, we investigate the relative contribution to the photoacoustic 
signal from the gold nanosphere and from its water environment, as a function of the pulse 
duration Tp and sphere radius Rg. The considered radii are on the order of a few nanometers 
to tens of nanometers, and the pulse durations typically range from tens of picoseconds to 
tens of nanoseconds. For each pair of parameters {Tp,Rs), two different numerical simula¬ 
tions were performed. Simulation (SI) computed the photoacoustic wave generated from 
the whole system, i.e the gold nanosphere and its water environment; simulation (S2) was 
identical except that the thermal expansion coefficient of gold was set to zero. The results 
from (S2) therefore only takes into account sound generation from water. Fig|2] shows a 
plot of the waveforms from simulations (SI) and (S2) for a 20-nm radius gold nanosphere 
for three different values of Tp (10 ps, 500 ps and 5 ns). It is clear from Figj2] that the 
predominant origin of the photoacoustic generation highly depends on the pulse duration: 
’’short” pulses mostly excite acoustic waves in the nanosphere, which are then radiated into 
the water, whereas photoacoustic waves with ’’long” pulses originate mostly from the liquid 
around the nanosphere. Both regimes have been studied experimentally. Indeed, various 
investigations have been conducted on acoustic vibration of gold nanoparticles in the short 
pulse regime (fs or ps excitation), see for instance^. In the nanosecond regime, Chen et 
ah have experimentally demonstrated that the photoacoustic signals originate from the 
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environment rather than the nanosphere itself^I. Their demonstration was based on the 
fact that the photoacoustic signal amplitude followed the properties of the temperature- 
dependence of the thermal expansion coefficient of the liquid around the particle. FigJ2] 
illustrates that our model and simulations encompass these different regimes, from the ex¬ 
citation of vibration modes in the sphere by short pulses (although no shorter than a few 
picosecond as a requirement of our thermal model) to photoacoustic generation directly in 
the surrounding liquid. It is therefore adapted to model a variety of different phenomena. 
In addition, our model can provide a quantitative assessment of the relative contribution to 



FIG. 3: Relative contribution from water to the overall photoacoustic energy (from water 
and gold), as a function of the laser pulse duration Tp and the nanosphere radius Rg. For a 
pulse duration of Tp = 5 ns, more than typically 90 % of the energy is emitted from water. 


the generated photoacoustic wave as a function of pulse duration Tp-. to do so, the energy e 
of the emitted photoacoustic wave (dehned as e = t)dt, r > Rg, independent 

of r) was computed for simulations (SI) and (S2). The relative contribution from water was 
dehned as r] = The values of p as a function of pulse duration and sphere diameter, 

plotted on FigjSl show that both the sphere diameter and the pulse duration affect the rel¬ 
ative contribution from water. However, for pulse durations larger than a few nanoseconds, 
most of the photoacoustic energy comes from the surrounding water, in agreement with 
the experimental results irt^^. For sphere diameters up to 40 nm, more than 90% of the 
photoacoustic energy is generated in water. 


From this point and throughout the rest of paper, we focus our interest on the nanosecond 
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regime, for which the photoacoustic emission from the nanosphere is negligible compared 
to that of water around it. Within this context, the following two sections investigate and 
quantify the typical dimension of the water layer that generates the photoacoustic wave, 
and compare the photoacoustic amplitude predicted for the nanospheres to those predicted 
from the point-absorber modeU^. 

2. Typical thickness of the generating water layer 



FIG. 4: Typical thickness of the water layer that emits the photoacoustic wave as a 
function of the pulse duration Tp, for different values Rs of the gold nanosphere radius. 

To quantify the size of the water layer that contributes to the photoacoustic generation, 
the following approach was implemented. The photoacoustic source term in EqJTO] may be 
straightforwardly turned off in the simulations by forcing to zero at any desired locations. 
Several simulations were therefore run by limiting the extent of the photoacoustic source 
term to distances r G [Rg] Rg + pgource]-, with pgource varied from 0 to -|-cxo. In practice, 
Psource was Varied up to a maximum value large enough so the photoacoustic signal did not 
differ signihcantly from its asymptotic value, corresponding to the case where all source 
points in water are active. The extent of the generating layer in water was then dehned 
by the value piayer = Psource for which the amplitude of the photoacoustic signal reached 
80 % of the amplitude of the asymptotic signal. This procedure was reiterated for different 
values of the laser pulse duration and the nanosphere radius, to compute the values of 
Piayeri^p, Rs) plotted ou FigJU FigJU shows that the size of the contributing layer is in 
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first approximation independent of the size of the sphere, and that it scales with the pulse 
duration approximately as piayer{Tp) ~ This scaling law suggests that the extent of 

the generating layer is dictated mostly by the diffusion of heat in water, regardless of the 
nanosphere diameter. As a consequence, each gold nanosphere may be considered as a 
nanometric absorber which thermally probes its environment within a spatial range driven 
by the laser pulse duration (longer than nanosecond). As an order of magnitude, a pulse 
duration = 5 ns yields player ~ 30 nm. 


3. Comparison with the photoacoustic point-absorber model 


It was shown above in Section IIII A II that for a nanometric sphere illuminated with a 
nanosecond pulse, the photoacoustic wave is mostly generated by the liquid surrounding the 
particle. Within this regime, the analytical model proposed by Calasso et al— for point- 
absorbers in the linear regime is therefore expected to predict reasonably well the amplitude 
and shape of photoacoustic waves generated by gold nanospheres. The objective of this 
section is to quantify the accuracy of this theoretical model by comparing its predictions to 
our simulations for hnite-size absorbers. This comparison will be further developed in the 
next results section for the nonlinear regime. From the analytical expression given by Eq. |5l 
the photoacoustic energy emitted by a point absorber is given by 

2 


e - f 

fcpomt,linear - ^abs ^3 ^3 J 


^+00 

[d/l 

/ —00 

du 


du (13) 

'^p 'p 

For gold nanospheres of different sizes, for a hxed pulse duration Tp = 5 ns, we com¬ 
pared the emitted photoacoustic energy predicted by Eq. [T3] to that predicted for gold 
nanosphere in the linear regime, for equivalent absorption cross-sections. Fig. |5] shows a 
plot of a = linear ^ functlon of the nanosphere radius Rs- The predictions from the 

^point,linear 

point absorber model and from our model turn out to be identical for vanishingly small 
diameters, as expected. Incidentally, this provides a validation of our numerical simulations 
in the linear regime. For hnite sizes, the effect of the gold nanosphere, both as an acoustic 
scatterer and as a photoacoustic source (via /3^), is to marginally decrease the emitted acous¬ 
tic energy compared to a point absorber of identical absorption cross-section. The effect is 
small, as expected from the very small ratio of the nanosphere diameter (typically tens of 
nm) to the acoustic wavelength (typically 20 pm for a 5-ns pulse). Therefore, in the nanosec- 
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FIG. 5: Ratio of the photoacoustic energy emitted from a gold nanosphere to that emitted 
from a point-absorber of identical absorption cross-section, as a function of the sphere 

dimensions, for a pulse duration Tp = 5ns. 


ond pulse regime, the point-absorber model proposed by Calasso et ai— provides accurate 
quantitative predictions for the emission of photoacoustic waves by a gold nanosphere, with 
overestimation of the acoustic energy less than 10 % for sphere radii up to 30 nm. 


B. nonlinear regime 

In this section, we quantitatively investigate for a gold nanosphere the impact of the 
temperature dependence of the thermal expansion coefficient, that leads to the so called 
thermal photoacoustic nonlinearity. In a preamble section, we first discuss the consequences 
that can be derived from the analytical expression provided by Calasso et ai— within 
the point-absorber model. We then report quantitative predictions obtained for a gold 
nanosphere in the nonlinear regime. 

1. Predictions from the analytic point-absorber model 

a. Existence of a critical absorbed energy. Eq. [7]shows that for fixed physical constants, 
the relative contribution of the nonlinear term to the photoacoustic pressure wave only 
depends on the absorbed energy Eats and the pulse duration Tp. For a fixed pulse duration 
Tp= 5 ns, Figj6] illustrates the change from the linear regime to the nonlinear regime as a 
function of the absorbed energy. We define the critical energy Ec as the value of absorbed 
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FIG. 6: Peak photoacoustic amplitude (measured at r = 1 mm) as a function of absorbed 
energy, computed from the point-absorber model by Calasso et a/.—, for an equilibrium 
temperature Tq = 20 °C and a pulse duration Tp = 5 ns. 



FIG. 7: Gritical energy Ec as a function of the equilibrium temperature Tq 


energy for which the peak amplitudes of the nonlinear and the linear contributions in Eq. [7] 
are identical. With this dehnition, the nonlinear contribution is thus signihcant for Eabs ^ 
Ec, and becomes predominant for Eabs ^ Ec- From Eq. [71 Ec is given by 


( 14 ) 
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The numerical prefactor can be computed numerically from the function h given by Eq. flA16j) 
in Appendix IA31 and yields 

Ec ~ 20.2 X PoCp§^(To)y^^ (15) 

Pi 

As indicated on Fig. El E,. = 31 fJ for Tp= 5 ns and Tq = 20 °C. By coupling Eq. IT^with 
the evolution of /9(T) in water shown in Fig. (H the temperature-dependence of the critical 
energy was computed and plotted on Fig. [3 As /9 vanishes around 4°C^ so does Ec, i.e. 
only nonlinear generation is predicted for any absorbed energy at that temperature. More 
importantly, the curve in Fig. [7] indicates that the critical energy is highly sensitive to the 
equilibrium temperature in a narrow range of a few degrees around 4°C. On the contrary, 
the critical energy is only weakly dependent on temperature at physiological temperatures 
(variation less than an order of magnitude for several tens of degrees). 

Interestingly, the critical energy given by the point-absorber model scales as the volume 
around the absorber over which heat diffuses during Tp. The criteria given by equation (3) 
by Egerev and Oraevsky— for a spherical particle of radius Rg can be restated as a critical 
energy with a form similar to that of the point-absorber model: 

Ec cx poCp^{To)Rl (16) 

Pi 

However, the volume term found by Egerev and Oraevsky— is the volume of the sphere, 
whereas it is given from Calasso et al.-^ by the volume of heat diffusion in water during 
the illumination pulse. These two predictions are incompatible: the critical energy given by 
Eq. [T6]goes to zero for vanishingly small spheres {Rg —)■ 0), whereas the point-absorber model 
predicts a hnite critical energy. Nevertheless, the two different expressions of the critical 
energy have in common to provide as a strong physical insight that the critical energy scales 
as some volume which origin remains to be determined. This will be further discussed based 
on our results for gold nanospheres in section IIIIB4I 

b. Frequency considerations. The difference in temporal shapes for the linear and the 
nonlinear contributions was pointed out by Calasso et al.—. Here, we investigate the con¬ 
sequences of these different temporal shapes (bipolar for the linear term, tripolar for the 
nonlinear term) in the frequency domain. For Tp = 5 ns. Fig. [S] shows the frequency spec¬ 
tra of the photoacoustic signal as a function of the absorbed energy, with each spectrum 
normalized by the absorbed energy. Fig. [8] shows that the spectrum amplitude varies non- 
linearly with the absorbed energy as expected from Eq. [7] in the time domain. Importantly, 


22 











0 50 100 150 200 250 300 350 400 

Frequency (MHz) 

FIG. 8: Photoacoustic amplitude spectra as a function of the absorbed energy, for 
Tp = 5 ns at To = 20 °C. Each spectrum was normalized by the corresponding absorbed 
energy Eats in order to illustrate the nonlinear dependence on Eats- 

it also indicates that the frequency content is shifted towards high frequencies in the non¬ 
linear regime. Therefore, the photoacoustic nonlinearity predicted by Calasso et al.^ only 
manifests itself for sufficiently high frequency. This is coherent with the fact that the critical 
energy decreases with decreasing pulse duration (increasing centre frequency), as indicated 
by Eq. [151 As a major consequence from the experimental point of view. Fig. |8] shows 
that even when the nonlinearity is predominant when considered over the full bandwidth, 
it remains minor for frequencies below 10 MHz even for Eats as high as 10 x E^. 

2. Temperature rise in a gold nanosphere 

It was shown in Section IIII A 31 in the limits of the linear model that the photoacoustic 
wave generated by a gold nanosphere is very close to that generated by a point absorber of 
identical absorption cross-section, and therefore depends only on the absorbed energy Eabs, 
with no dependency on thermal properties. In the case of thermal nonlinearity (caused by 
the temperature-dependence of the thermal expansion coefficient), it is the temperature held 
that drives the effective value of (3{T). As a consequence, because the temperature helds are 
different for a point absorber and a gold nanosphere with the same absorption cross-section, 
one expects the thermal photoacoustic nonlinearity to be dependent on the size of the ab¬ 
sorber. As opposed to the point-absorber model, all the values of the temperature held are 
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FIG. 9: Temperature rise inside a 20-nm radius gold nanosphere illuminated with a 
fluence <I)o = 5mJ.cm“^ {Eats ~ 165 fJ) and Tp = 5 ns. 

hnite when hnite-size absorbers such as gold nanospheres are considered. Moreover, for a 
given incident fluence, the temperature rise in a nanosphere is highly dependent on its size, 
both through the size-dependence of the thermal diffusion (the thermal Green’s function 
given in Appendix |B] depends on the particle radius) and the absorption cross-section (see 
table ttH). Fig. [9] illustrates the temperature rise in a 40-nm diameter gold nanosphere illumi¬ 
nated with a fluence $o = 5mJ.cm“^ {Eabs ~ 165 fJ) and Tp = 5 ns. This plot shows that 
for a fluence value typical of those used for biomedical applications, the temperature rise in 
a 40-nm diameter gold nanosphere is signihcantly larger than the equilibrium temperature. 
In this case, the photoacoustic nonlinearity is likely to become signihcant, as demonstrated 
further below. For a given illumination fluence. Fig. [10] shows that the peak temperature 
rise depends on the sphere radius and turns out to be maximum for a radius value around 35 
nm, reflecting the dependence on radius via both thermal diffusion and the size-dependence 
of the absorption cross-section. 


3. Gold nanosphere vs. point absorber 

Simulations were run to predict the photoacoustic signals generated from gold nanospheres 
of various diameters illuminated with 5-ns pulses of various fluences, in order to compare the 
photoacoustic amplitude to that predicted from the point-absorber model^ with matched 
absorption cross-sections. Fig. ITT] shows a plot of the ratio of the photoacoustic peak am- 
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FIG. 10: Peak temperature rise as a function of gold nanosphere radius, for a fluence 
$0 = 5mJ.cm"2 [Eabs 165 fJ) and = 5 ns. 



FIG. 11: Photoacoustic amplitude from a gold nanosphere, normalized to that of the 
point-absorber model, as a function of nanosphere diameter. Parameters: Tq = 20 °G, 

Tp = 5 ns 

plitude predicted for the nanosphere to that predicted for the point-absorber as a function 
of the absorbed energy. Figure [TT] confirms that in the linear regime, i.e for low enough 
absorbed energy (or equivalently Eabs "C Ec), a gold nanosphere may be considered as a 
point-absorber, i.e. the predictions from the corresponding point-absorber model are accu¬ 
rate. However, in the nonlinear regime, the point-absorber model significantly overestimates 
the amplitude of the photoacoustic signals. In other words, the critical energy for a gold 
nanosphere is significantly higher than that predicted by the point-absorber model. For the 
gold nanosphere, in order to dehne the critical energy in agreement with the dehnition for 
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the point absorber (see sec lIII B 1 al) . the nonlinear contribution was dehned by the difference 
between the total signal predicted in the nonlinear regime and the signal predicted in the 
linear regime only (by keeping /3 constant). As an illustration, the photoacoustic amplitude 



FIG. 12: Photoacoustic peak amplitude for a 20-nm radius gold nanosphere and a 
point-absorber of identical absorption cross-section. Tq = 20 °C and Tp = 5 ns 

predicted for a 20-nm radius gold nanosphere at equilibrium temperature Tq = 20 °C, illu¬ 
minated with a 5-ns pulse of fluence $0 = 5mJ.cm“^ {Eabs ~ 165 fJ), is about three times 
lower than that predicted by the point-absorber model, as illustrated in Fig. [121 ^md the 
critical energy is 216 fJ whereas the point-absorber model predicts a value of 31 fJ. On the 
other hand, the frequency features discussed in the previous section for the point-absorber 
model remain strictly identical for a gold nanosphere, i.e. the nonlinearity is observed 
only at high enough frequencies. Fig. |H] for the point-absorber remains valid for a gold 
nanosphere, provided that the appropriate value of critical energy is used. In conclusion to 
this comparison between our predictions for gold nanospheres and predictions for the point- 
absorber model, the hnite size of the gold nanosphere must be taken into account to obtain 
accurate quantitative predictions regarding the occurence of thermal-based photoacoustic 
nonlinearities. 

4- Optimal size for nonlinear generation 

In order to quantitatively predict the occurence of nonlinear photoacoustic generation 
as a function of parameters that are controllable experimentally, simulations were run for 
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TABLE III: Values of the critical energy Ec, critical fluence <I>c and maximum temperature 
rise ATmax for gold nanospheres as a function of radius Rg and equilibrium temperature 

To, with Tp = 5 ns. 


To{°C)\Rg (nm) 


10 

20 

30 

40 

50 

60 


A«^(fJ) 

9.4 

19 

38 

66 

no 

170 

5 

4>^®(mJ/cm^) 

2.5 

0.59 

0.35 

0.36 

0.51 

0.74 


ATmax{°C) 

17.4 

13.7 

13.8 

14.5 

15.5 

17.1 


A«^(fJ) 

52 

93 

150 

230 

340 

490 

10 

4>^''(mJ/cm^) 

14 

2.8 

1.4 

1.3 

1.6 

2.1 


ATmax{°C) 

96 

65 

55 

51 

49 

49 


A«^(fJ) 

120 

220 

350 

520 

740 

1040 

20 

4>^®(mJ/cm^) 

33 

6.5 

3.2 

2.9 

3.6 

4.5 


ATmax{°C) 

230 

150 

126 

114 

107 

104 


Af^(fJ) 

230 

400 

630 

940 

1350 

1850 

40 

4>^'*(mJ/cm^) 

62 

12 

5.8 

5.1 

6.4 

8.1 


ATmax{°C) 

420 

280 

230 

210 

190 

190 


gold nanospheres with different radii and equilibrium temperatures. Table. IIIII reports the 
corresponding results as the value of the critical energy Ec, critical fluence <I>c = Ec/da 
and the peak temperature in the sphere. In particular, these results show that for any 
equilibrium temperature, there is an optimal sphere radius around 40 nm for which the 
critical fluence is minimal (as a function of size). In other words, our model predicts that 
at a given illumination fluence, the thermal-based photoacoustic nonlinearity is maximised 
for a sphere radius around 40 nm. We note that the value of the gold nanosphere radius 
that maximizes the photoacoustic nonlinearity (for a given illumination fluence) is close to 
the one that maximizes the peak temperature rise in the nanosphere (see Fig. fTOj) : because 
our model assumes the continuity of temperature across the gold/water interface, the peak 
temperature in water is also maximized for a radius of 40 nm, which is expected to maximize 
thermal nonlinearity. The values of critical energy given in Table. HTTI can be further analyzed 
as a function of the particle size, for a given equilibrium temperature, in order to assess the 
relevance of the scaling law predicted by the point-absorber model {Ec independent of Rg, see 


27 















FIG. 13: Critical energy as a function of nanosphere radius (logarithmic axis). 

Eq. US]) and the model by Egerev and Oraevsky— [Ec oc see Eq. [16]). For Tq = 20 
Fig Us] shows a plot of the values of the critical energy predicted by our model, by the 
point-absorber model, and for Ec oc and Ec oc This plot is compatible with the 
hypothesis that the critical energy tends towards the constant value of the point-absorber 
for small spheres, and scales as the surface for large spheres. For large spheres, it is 
indeed expected that the volume of heat diffusion scales as the surface of the sphere times 
the radial distance traveled by heat (only dependent on the pulse duration and the diffusivity 
in water). Beyond providing quantitative values for the critical energy (or equivalently the 
critical fluence via the absorption cross-section), our model therefore also provides some 
physical insight into the origin of the critical energy. 


5. Influence of the equilibrium temperature 

In this section, we study the influence of the equilibrium temperature on the photoacous¬ 
tic amplitude, for a fixed illumination fluence. Fig. [TT] shows the results obtained in the 
linear regime (lower black curve) and in the nonlinear regime (upper blue curve). In the lin¬ 
ear regime, the photoacoustic amplitude reflects the temperature dependence of the thermal 
expansion coefficient /9o(To), as expected from Eq. |2l In the nonlinear regime, at a higher 
illumination fluence, the temperature-dependence of the signal is much less signihcant, and 
essentially reflects the effective value of the thermal expansion coefficient, different from the 
value at equilibrium because of the signihcant temperature rise around the particle. As a 








FIG. 14: Peak photoacoustic amplitude emitted by a 40-nm diameter gold nanosphere, as 
a function of equilibrium temperature, for two different fluence 4). 4> = 0.1 mJ/cm^ 
corresponds to the linear regime, whereas signihcant nonlinearities occur at 

$ = 10 mJ/cm^. 


consequence, the evolution of the peak photoacoustic amplitude as a function of equilibrium 
temperature reflects temperature rises at the nanoscale. In particular, the strong signal 
observed at Tq = 4 °C for which /3o = 0 is purely nonlinear in nature as the linear contri¬ 
bution vanishes at this temperature. This has hrst been observed experimentally by Hunter 
et al.— in protons experiments and has been more recently reported with gold nanoparticles 
by Fukasawa et al.— and Simandoux et al.—. Recent experimental measurements of the 
photoacoustic amplitudes as a function of equilibrium temperature with 40-nm diameter 
gold nanospheres by Simandoux et al.— are in good agreement with the predictions given 
in Fig. [m 

An important consequence of the results shown in Fig. [Tl]is that when thermal nonlin¬ 
earity takes place, the amplitude of photoacoustic waves cannot provide a measurement of 
the equilibrium temperature through the knowledge or calibration of /9o(R))- While pho¬ 
toacoustic measurements can provide a measurement of the equilibrium temperature in the 
linear regime, as was originally proposed by Larina et al.— and further investigated in sub¬ 
sequent works (see—-— for instance), this approach is expected to not work any longer with 
nanoparticles illuminated with high enough fluences. 
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C. Discussion 


The photoacoustic amplitude predicted by the point absorber model (see Eq. E]) is inde¬ 
pendent of thermal diffusion in water, and only depends on the absorbed energy. Therefore 
our results suggest that when gold nanospheres behave as point-absorbers in the linear 
regime, the photoacoustic amplitude is insensitive to heat transfer from gold to water, 
although the temperature field does depend on the thermal diffusivity of water (see Ap¬ 
pendix |B|). On the contrary, nonlinear photoacoustic signals are strongly dependent on the 
thermal diffusivity of water (as predicted by the point-absorber model in water by Eq. [7]). 
This suggest that the occurence of thermal nonlinearity is required to probe thermal diffusion 
properties at the spatial and temporal nanoscale. In the model described here, the tempera¬ 
ture was assumed to be continuous across the gold/water interface, whereas it is known that 
an interfacial thermal resistivity exists and may have significant impact at the nanoscale^^— . 
To check the importance of this effect on our predictions, we solved the thermal problem 
by taking into account an interfacial thermal resistance at the gold-water interface. The 
temperature held for a delta impulse illumination (Green’s function) was found to strongly 
depend on the thermal resistance, in agreement with earlier work—. However, when con¬ 
volved with a 5-ns pulse illumination, the temperature held for a interfacial resistance of 
10“® m^.K.W”^, typical of the gold-water interface, was nearly identical to the case with no 
thermal resistance (less than a few % relative diherence). The results presented in this work 
with no interfacial thermal resistance therefore remains valid in the few-nanosecond regime. 
On the other hand, the presence of coating, such as silica coatings or PEG (poly-ethylene 
glycol) coatings, may have a signihcant ehect, and will be investigated in a future work. 

Importantly, our work describes the photoacoustic generation by a single gold nanosphere. 
When collections of nanoparticles such as encountered in suspension are considered, the 
emitted photoacoustic waves arise from the sum of each contribution from individual par¬ 
ticles. In this case, the frequency content of the resulting wave is dictated by the spatial 
distribution of the collection of nanoparticles, which acts as a low-pass filter. Because the 
thermal nonlinearity was shown to be significant only at high frequencies, it is expected 
that it may not be observable at low frequencies. In a recent experimental study^ with a 
detection frequency of 20 MHz and 40-nm diameter gold nanospheres suspended in 100 /xm 
diameter tube, a linear relationship was observed between the peak-to-peak amplitude and 
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the fluence, up to a fluence of <h = 7 mJ/cm^, yet close to the corresponding critical flu- 
ence predicted here <hc = 6.5 mJ/cm^. However, as discussed above, the critical fluences 
predicted here are only valid for a single nanoparticle, and the prediction for ensembles of 
nanoparticles, beyond the scope of this work, requires taking into account the nanoparticles 
spatial distribution. 

It also has to be kept in mind that our model assumes that there is no limitation on 
the peak temperature in the gold nanosphere and in water. In practice, the model becomes 
of course irrelevant if the predicted temperatures would lead to phase transitions in either 
gold or water. In addition, the thermal model used in our study is strictly restricted to 
the long-pulse illumination regime (pulses typically longer than a few picosecond), and the 
corresponding modeling cannot be applied to the sub-picosecond illumination regime which 
usually requires either a two-temperature model^ or nonthermal distributions^ii^^ for the 
electrons and phonons in nanoparticles. 


IV. SUMMARY AND CONCLUSION 

In summary, we theoretically investigated the photoacoustic generation by a gold 
nanosphere in water in the thermoelastic regime. Photoacoustic signals were predicted 
numerically based on the successive resolution of a thermal diffusion problem and a ther¬ 
moelastic problem, taking into account the hnite size of the gold nanosphere, thermoelastic 
and elastic properties of both water and gold, and the temperature dependence of the 
thermal expansion coefficient of water. For sufficiently high illumination fluences, this 
temperature-dependence yields a nonlinear relationship between the photoacoustic ampli¬ 
tude and the fluence. For nanosecond pulses in the linear regime, we showed that more 
than 90 % of the emitted photoacoustic energy is generated in water, and the thickness of 
the generating layer around the particle scales close to the square root of the pulse dura¬ 
tion. In the linear regime, we showed that the point-absorber model introduced by Calasso 
et al.^ accurately predicts the amplitude of the photoacoustic waves generated by gold 
nanospheres for diameters up to several tens of nanometers. However, the point-absorber 
model signihcantly overestimates the amplitude of photoacoustic waves generated by gold 
nanospheres in the nonlinear regime. Our model for hnite-size particles provided quanti¬ 
tative estimates of the critical energy, dehned as the absorbed energy required such that 
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the nonlinear contribution is equal to that of the linear contribution. Our results suggest 
that the critical energy scales as the volume of water over which heat diffuses during the 
illumination pulse. A frequency analysis of the nonlinear signals indicated that the ther¬ 
mal nonlinearity from a gold nanosphere is more pronounced at high frequencies dictated 
by the pulse duration. Finally, we show that the relationship between the photoacoustic 
amplitude and the equilibrium temperature at sufficiently high fluence reflects the thermal 
diffusion at the nanoscale around the gold nanosphere. Although our model was limited 
to the case of a bare nanoparticle, the approach used in this work can be extended to the 
more general case of a coated nanoparticle, provided that the corresponding temperature 
held can be predicted, either by analytical or numerical means (such as given by Baffou 
and Rigneault— ). To limit the computational costs, and for sake of simplicity, this work 
was limited to nanospheres with central symmetry. However, the same methodology could 
apply in two dimensions for axisymmetric nanoparticles such a gold nanorods at the cost of 
more intensive computations. 
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Appendix A: Nonlinear thermoelastic regime 


1. Derivation of Eq. [B] 

We consider the linear regime of acoustic propagation, i.e. small Mach number (^ 1) 

and small density changes (^^ 1). In particular, the following approximation holds: 

^ ^ + (v.V)X ~ The three fundamental equations of linear acoustics can be 
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written as 


^ = -Po div(v) 

(Al) 

dv 

(A2) 


ds 

pT— = kAT + Py 

(A3) 


where p{r,t) is the acoustic pressure, v(r,t) is the acoustic displacement velocity, s(r,t) is 
the specihc entropy and T{r,t) is the temperature. The thermal conduction coefficient k 
is assumed to be constant. Moreover, the differentials of the state functions p = p{p,T) et 
s = s{T,p) can be written as^^ 


7 


Sp=MSp-flST] 


fc = ^|iT+'^ 


7/3 


where the thermodynamical coefficients are dehned as 


-6p] 


(A4) 

(AS) 


c? = 


Cp = T 


dp 

dp 

ds 

&T 


c„ = T — 


(S), 


7 = — 

Cy 

with the isentropic sound velocity, Cp the specihc heat capacity at constant pressure, Cy the 
specihc heat capacity at constant volume and /3 the thermal expansion coefficient. Eqs. 
and IASI may be written as 


^ _ R^] 

dt cl dt dt 

ds Cp,dT '^ — Idp 

dt T dt ^ '-^13 dt 


(A6) 

(A7) 


where all the thermodynamical coefficients may vary with properties such as temperature. In 
the following, we consider that only the variations of /3(T) with temperature are signihcant, 
and assume all other thermodynamical properties to be constant. By coupling the hrst two 
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fundamental equations lAll and IA21 one gets — Ap = 0, which can be transformed into 


the following equation by use of Eq. 

7 d'^p 


— Ap = 


I3{T)T 


ds 

m 


(A 8 ) 


cl dt^ Cp dt 

By use of Eq. flA7p into the evolution equation 
equation for the temperature held: 

dT 

PO^V o, KAT -|- I Q rN, 

at 7/3 at 

Under the assumption, valid for liquids, that 7 ~ 1, one eventually obtains Eqs. (II]) and d?]) 
given in section III Al 


for the entropy, one obtains the following 
^ ^ (A9) 


dT 

POCy— = KAT + Py 
at 


1 d'^p 


d 


Ap = po—J/3{T)— j 


dT\ 


(AlO) 

(All) 


cl dP dt 

Eq. (lAlip may also obtained from Eqs. (IT 0 |) . by considering a liquid as a special type of 
solid for which p = 0 . 


2. Consequences from the form of Eq. dTj) / (|A11() 


The photoacoustic wave equation Eq. d?]) / (lAlip used in this work and derived above, 
differs from the equation originally introduced by Burmistrova et al.— which is given only 
for small temperature variations around Tq and reads: 


= (A12) 

Eq. (1A12P corresponds to that given by Gusev and Karabutov-, further used by Calasso 
et al— for the analytic derivation corresponding to the point-absorber model. Although the 
two forms source term are very different in these two equations, it turns out as demonstrated 
below that the consequences on the solution are relatively minor when one consider only 
the hrst-order dependence of (3{T) with temperature, as /5(T) = (3o + (3i{T — Tq) (with 
/So = /S(So) and /Si = ^(Tq)). By using this hrst-order expansion of /S(T) around Tq, our 
equation Eq. (CP becomes 


1 d‘^p 
cl dP 


- Ap = po 


{/3o + (3i (T 


d^T 



(A13) 
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whereas the original equation used by Gusev and Karabutov-, Calasso et a/.—, Burmistrova 
et ai— becomes 


1 d'^p 
cl dt"^ 


- Ap = po 


{/3o + 2/?i(T 


d^T 

To)) ^ +2ft 



(AM) 


It therefore turns out that the only remaining difference is a factor 2 before /3i. In conse¬ 
quence, all the analytical calculations based on Eq. flAl2ll . as given by Gusev and Karabu¬ 
tov-, Galasso et al.—, Burmistrova et al.—, may be used straightforwardly by changing (3i 
to /3i/2. 


3. Modified expression for the point-absorber model 

Three modihcations are needed to go from the original equation (25) given by Galasso 
et ai— to our Eq. 

1. Galasso et al— erroneously used a term in their Eq. (19), instead of as 

straightforward from their Eq. (6). As can be derived from simple calculations, the 
only consequence of this minor error is that equation (25) by Galasso et al.— has to 
be multiplied by a factor 2^/^. 

2. In Galasso et al.—, the derivation of the nonlinear contribution to the photoa¬ 

coustic signal by a point-absorber was based on Eq. (I A 121) (or stricly speaking 
its equivalent formulated as a displacement potential), under the assumption that 
/S(T) = + (3i{T — To). From the discussion above, we therefore had to change 

the expression given in Galasso et al— by replacing /3i to /5i/2 to take into account 
our equation [71 (N.B. Galasso et ai— used the notation {(31,(32) whereas we used 
{(30,(31)). In other words, Eq (25) in Galasso et al— simply has to be multiplied by a 
factor 1/2 to take into account Eq. ([7]), but all the analytical derivations remain valid 
except for this numerical prefactor. 

Overall, to take into account the two prefactors above, Eq (25) in Galasso et al.— has 
to be multiplied by a factor 1/2 x 2^/^ = 2^/^. 

3. The pulse duration 6 used by Galasso et ai— is related to the pulse duration Tp 
dehned here by Tp = 2^Jla{2) 9. Moreover, we use the dimensionless retarded time 


35 











f = {t — r/Cs)/Tp whereas Calasso et al— used r = {t — r/cs)/6. This change of 
variable further adds a numerical prefactor to the expression of the nonlinear term. 


Taking into account the three modihcations above, the nonlinear contribution from the 
point-absorber model can be written as 




47rr 


1 ( t- - 

h f= -^ 


Tn 


where h{f) is given by 


h{f) = 


[ln(2)]3/4 d^ 
Stt^ dr^ 


erf(^/ y/2) ^-2fg/2-2^1n(2)). 

^3/2 


de 


(A15) 


(A16) 


Appendix B: Green’s function of the thermal problem 


The following expression of the Green’s function of the thermal problem may be found 
in Egerev et al.— and written as 


where: 


Gth(r,t) = 


i?o 


«! — 02 r 


+ ai exp ( ai(^-1) -f- a\— ) erfc * 


02 exp ( ot 2 {— -1) + ot 2 — ) erfc 


R, 


'^th 



Rl 


V 


Tth = - , ai = - 1 ± W 1 - - , V = ^- An An 

X 2 2 V V vJ Po Cp 


Po S 


(Bl) 


Appendix C: FDTD algorithm for the thermoelastic problem 


Following the approach introduced for elastodynamics by Virieux— , Eqs. [TO] were dis¬ 
cretized on staggered grids (see Fig. ITB]) with the following centered finite-difference approx¬ 
imation: 

(Cl) 


da Aa 

where a may be either r or t. Moreover, whenever needed at locations where they were 
undehned (either held variables or material properties), variables were approximated by 
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their arithmetic mean with /(a*) = ^ )^/K — 2 _}_^ material properties were dehned 

at locations kAr. The Eqs. [lO] accordingly read 


Vr{k + |,m + i) — Vr{k + — i) 

At 

arr{k + 1, m) — arr{k, m) 


po{k + 1) + Po{k) 


X 


Ar{k + i) 


Ar 


X 


arrik + l,m) + arr{k,m) aee{k + l,m) + aee{k,m) 


(C2a) 


arr{k, m + 1) — arr{k, m) 


At 




, /, 1 In fi 1 In 

k + -\ Vr[k + -,m+-)- \k--\ Vr[k- -,m + -) 


+ 


2p{k) 


Ar 


’ 1 In 1 In' 

Vr{k + -,rn+ -) -Vr{k- -,m+ -) 


Kk) + ^^k) 




(C2b) 


aee{k, m + 1) - aee{k, m) 


At 




, /, 1 In a 1 In 

k + -j Vr{k +-,m+- \k --j Vr{k- -,m + -) 


+ 


2p{k) 


1 1 
Ar 2 


Vr{k + \,m + \) ^ Vr{k - \,m + \y 


k + \ 


+ 


k — h 


^{k) + -p{k) 




(C2c) 


where we used (^ + f) At the central point, Eqs. IC2b and IC2b cannot be 


used because of the singularity in 1/r. At r = 0, we used ^ ~ r^ and ^^(0) = 0 to obtain 
the following form 


^^(0, t) = (3A(0) + 2^(0)) A(0, i) - (a( 0) + /3(r)^(0, t) (C3a) 




(C3b) 


By using |^(0,f) the following discretized equations are hnally obtained for 


Vr{lAr) 
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At 


r = 0 


j 

1 n/ 
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FIG. 15: Spatio-temporal mesh used in FDTD simulations. The stress held arr (blue 
points) and the radial velocity displacement (red arrows) v are discretized over staggered 
grids, both in time and spaced. The grids are staggered such that pressure variables are 
dehned on the central point r = 0. Values of aee and T are dehned at the same positions 

3jS • 


the central point: 

(Jrr(0, m -|- 1) — arr(0, m) 

At 


= [3A(0) + 2/i(0)] m + 


-^(0) + 2^(0) 


/3(T(0,m + l/2)) 


T(0, m -1-1) — T(0, m) 


At 


ae 0 (O,m) = arriS^.m) 

The temporal step was to the spatial step by the following stability condition: 

Ar 

At = 0.99 X 


(C4a) 

(C4b) 


Cgold 

with Cgoid = 3240 m/s the highest speed of sound involved in the problem. 
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